Event-by-event fluctuations of azimuthal particle anisotropy 
in Au+Au collisions at ^/s^= 200 GeV 
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This paper presents the first measurement of event-by-event fluctuations of the elliptic flow 
parameter V2 in Au-|-Au collisions at ys„„ = 200 GeV as a function of collision centrality. The 
relative non-statistical fluctuations of the V2 parameter are found to be approximately 40%. The 
results, including contributions from event-by-event elliptic flow fluctuations and from azimuthal 
correlations that are unrelated to the reaction plane (non-flow correlations), establish an upper limit 
on the magnitude of underlying elliptic flow fluctuations. This limit is consistent with predictions 
based on spatial fluctuations of the participating nucleons in the initial nuclear overlap region. These 
results provide important constraints on models of the initial state and hydrodynamic evolution of 
relativistic heavy ion collisions. 
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Results from the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory suggest that 
a dense state of matter is formed in ultrarelativistic 
nucleus-nucleus collisions [IHl]. Studies of final state 
charged particle momentum distributions show that the 
produced matter undergoes a rapid collective expansion 
transverse to the direction of the colliding nuclei. In 
particular, for collisions at non-zero impact parameter, 
the expansion shows a significant anisotropy in the 
azimuthal angle, strongly correlated with the anisotropic 
shape of the initial nuclear overlap region. The dominant 
component of this anisotropic expansion is called "elliptic 
flow" and is commonly quantified by the second coeffi- 
cient, V2, of a Fourier decomposition of the azimuthal 
distribution of observed particles relative to the event- 
plane angle [5]. 

Elliptic flow has been studied extensively in Au+Au 
collisions at RHIC as a function of pseudorapidity, 
centrality, transverse momentum and center-of-mass en- 
ergy PHllISHH]. For Au+Au coUisions at RHIC energies, 
the observed dependence of the elliptic flow signal on 
centrality and transverse momentum is found to be 
in good agreement with calculations in hydrodynamic 
models [SI [H]. This is considered evidence for an early 
equilibration of the colliding system and a low viscosity 
of the matter produced in the early stage of the collision 



process [10]. In such calculations, for given conditions of 
the produced matter, the elliptic flow magnitude is found 
to be proportional to the eccentricity characterizing the 
transverse shape of the initial nuclear overlap region |llj . 

Measurements of elliptic flow in the smaller Cu+Cu 
system have shown surprisingly large values of elliptic 
flow, in particular for the most central collisions where 
the average eccentricity of the nuclear overlap region was 
expected to be small |12j . Experimental measurements 
of V2 can be affected by event-by-event fluctuations in 
the initial geometry [13] and it is possible to reconcile 
the results for Cu+Cu and Au+Au collisions if these 
fluctuations are properly accounted for [T^]. To this 
end, we have proposed a new deflnition of eccentric- 
ity, which does not make reference to the direction of 
the impact parameter vector, but rather characterizes 
the eccentricity through the event-by-event distribution 
of nucleon-nucleon interaction points obtained from a 
Glauber Monte-Carlo calculation [121 E]. This method 
of calculating the initial state anisotropy, which leads to 
flnite "participant eccentricity" values even for the most 
central events and has a large effect in the smaller Cu+Cu 
system, has been found to be crucial for understanding 
the comparison of Cu+Cu and Au+Au elliptic flow 
results 

Using the probabilistic distribution of interaction 



2 



points obtained from a Glauber calculation, performed 
on an event-by-event basis, leads to relative eccentricity 
fluctuations of (J^^^^^/{epaTt) «40% for Au-|-Au collisions 
at fixed number of participating nucleons (iVpart) [H]- 
Similar calculations taking into account stochastic initial 
state interaction points in a color glass condensate (CGC) 
model also yield large relative eccentricity fluctuations 
of '''epart/i^part) ~ 30% [U] . If V2 is proportioual to e, 
an event-by-event measurement of elliptic flow should 
therefore exhibit sizable fluctuations in V2, even at fixed 

iVpart. 

An event-by-event measurement of the anisotropy in 
heavy ion collisions is expected to yield fluctuations from 
three sources: statistical fluctuations due to the flnite 
number of particles observed, elliptic flow fluctuations 
and other many-particle correlations. The statistical 
fluctuations in the observed V2 signal can be taken out 
with a study of the measurement response to the input 
V2 signal. Particle correlations other than flow (non-flow 
correlations) such as HBT, resonance decays and jets can 
resemble correlations due to elliptic flow and have various 
effects on different flow measurements. In particular, 
non-flow correlations can broaden the apparent V2 dis- 
tribution and enhance the observed V2 fluctuations. This 
letter presents the flrst measurement of event-by-event 
dynamic fluctuations in V2, which include contributions 
from elliptic flow fluctuations and non-flow correlations. 

The data shown here were taken with the PHOBOS 
detector at RHIC during the year 2004. The PHOBOS 
detector is composed primarily of silicon pad detectors 
for tracking, vertex reconstruction, and multiplicity mea- 
surements. Details of the setup and the layout of the 
silicon sensors can be found elsewhere [16]. The collision 
trigger, event selection and centrality determination are 
described in Ref. [T7]. The Monte Carlo (MC) simu- 
lations of the detector performance are based on the 
HIJING event generator [18] and the GEANT 3.21 [H] 
simulation package, folding in the signal response for 
scintillator counters and silicon sensors. 

The PHOBOS multiplicity array, composed of single 
layer silicon pad detectors, has a coverage of \ri\ < 5.4 
over almost the full azimuth. We parametrize the pseu- 
dorapidity dependence, f2(??), with a single parameter, 
V2 = i'2(0), and a triangular or trapezoidal shape, given 

by -r(^) = V. (1 - M), or (,) = {|:jl;i<:^,„>„ 

respectively. Both of these parameterizations provide a 
reasonable description of the measured pseudorapidity 
dependence of elliptic flow [5] . 

The event-by-event measurement method has been 
developed to use all the available information from the 
multiplicity array to measure the elliptic flow at zero 
rapidity, V2 , while allowing an efficient correction for the 
non-uniformities in the acceptance. Taking into account 
correlations only due to elliptic fiow, the probability of 
a particle with given pseudorapidity, 77, to be emitted 
in the azimuthal angle, 0, in an event with elliptic fiow 



magnitude, V2, and event-plane angle 0o is given by 

pi(j)\v2, 0o; ?7) = ^ {1 + 2v2{v) COS (2 [(j) ~ (^0])} • (1) 

The direction of the event-plane angle, (j^o, is expected to 
align with the reaction plane angle if the initial geometry 
of heavy ion collisions is defined by two smooth Wood- 
Saxon distributions or with the participant eccentricity 
axis if initial geometry fluctuations are indeed present. In 
this measurement, 4>o is determined from the distribution 
of final state particles without relying on any model 
about the initial geometry of the collision. 

The angular coordinates (77, (j)) of charged particles are 
measured using the location of the energy deposited in 
the silicon multiplicity detectors. After merging of sig- 
nals in neighboring pads in cases where a particle travels 
through more than a single pad, the deposited energy 
is corrected for the angle of incidence, assuming that 
the charged particle originated from the primary vertex. 
Noise and background hits are rejected by placing a 
lower threshold on this angle-corrected deposited energy. 
Depending on 77, merged hits with less than 50-60% of 
the energy loss expected for a minimum ionizing particle 
are rejected [2D]. Since the multiplicity array consists of 
single-layer silicon detectors, there is no px, charge or 
mass information available for the particles. All charged 
particles above a low-px cutoff of about 7 MeV/c at 
?7=3, and 35 MeV/c at 77=0 (the threshold below which 
a charged pion is stopped by the beryllium beam pipe) 
are included on equal footing. We define the probability 
density function (PDF) for a hit position (77, cf)) for an 
event with V2 and event-plane angle 4>o as 

P{<l>\v2,(l)o;v) = \ rP('^|v2, 00; v), (2) 

S(V2,(po;Vj 

where the normalization parameter s(y2, 4'o]''l) is calcu- 
lated in small bins of 77 such that the PDF folded with the 
acceptance is normalized to the same value for different 
values of V2 and (f>Q. The normalization parameter is 
given by 

rri+Ari 

s(w2>o,?7) = / ^(??','/')p('/'l^'2,0o;?7')d0d77', (3) 

where the acceptance function, ^(77, (jj) denotes the prob- 
ability of a particle moving in the 77, (jj direction to yield 
a reconstructable hit. 

For a single event, the likelihood function of V2 and (po 
is defined as L{v2,(t)Q) = 117=1 ^(0j|v2, '/'o; »7i), where the 
product is over all n hits in the detector. The likelihood 
function describes the probability of observing the hits in 
the event for the given values of the parameters V2 and 
4>Q. The parameters V2 and (f>Q are varied to maximize 
the likelihood function and estimate the observed values, 
V2^^ and (I)q^^, for each event. 

The response of the event-by-event measurement is 
non-linear and depends on the observed multiplicity n. 
Therefore, a detailed study of the response function is 
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required to extract the true V2 distribution from the 
measured V2^^ distribution. Let /(V2) be the true V2 
distribution for a set of events in a given centrahty bin, 
and givf^^) the corresponding observed distribution. The 
true and observed distributions are related by 

5(vf ^) = / i^(vf ^ V2, 71) /(V2) dv2 N{n) dn, (4) 



where N(n) is the multiphcity distribution of the given 
set of events and K{v2^^,V2,n) is the expected distri- 
bution of 1^2'^'' for events with fixed input flow V2, and 
constant observed muhiphcity n. 

The response function, K{v2^^ ,V2,n) is determined 
by performing the event-by-event analysis on modified 
HIJING events with fiow of fixed magnitude V2. The 
flow is introduced by redistributing the generated par- 
ticles in each event in the (f) direction according to 
the probability distribution given by Eq. [T] and the 
assumed pseudorapidity dependence of V2- For the two 
parameterizations of V2{'ri), triangular and trapezoidal, 
used in the event-by-event measurement, the correspond- 
ing response functions, K^" and K^^'^^p, are calculated. 
Fitting smooth functions through the observed response 
functions decreases bin-to-bin fluctuations and allows for 
interpolation in V2 and n. The response of a perfect 
detector can be determined as a function of event multi- 
plicity In practice, some empirical modiflcations 
to the ideal relation, accounting for detector effects, 
significantly improve fits to the response function, leading 
to 
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FIG. 1: (?;2) (top) and a^yn (bottom) versus A^part for 
Au-|-Au collisions at ys„„ = 200 GeV. Previously published 
event-plane V2 results for the same collision system are shown 
for comparison Boxes and gray bands show 90% C.L. 
systematic errors and the error bars represent 1-cr statistical 
errors. The results are for < r; < 1 for the track-based 
method and \ri\ < 1 for hit-based and event-by-event methods. 
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with v^""* ^ (An + B)v2 and a = C/^n + £>, and 
where /q is the modified Bessel function. The four 
parameters {A, B, C, D) are obtained by fits to observed 
K{vf'',V2,n) in the modified HIJING samples. 

Correcting for all known effects incorporated in our 
MC, we obtain the true event-by-event V2 distribution, 
/(V2), which includes contributions from elliptic flow 
fluctuations and non-flow correlations. We assume /(V2) 
to be a Gaussian in the range V2 > [Hj with two 
parameters, V2 and , denoting the mean and standard 
deviation in the given range. For given values of the 
parameters, it is possible to take the integral in Eq. |4] 
numerically to obtain the expected V2'"' distribution. 
Comparing the expected and observed distributions, the 
values of V2 and cr,, are found by a maximum-likelihood 
flt. Midrapidity {\irj\ < 1) results from the two parameteri- 
zations of W2(?7), triangular and trapezoidal, are averaged 

0.5(iiv2*" + V2*"''P) and 



> 12 
0.5(iicr,*" 



trap) ^j^g 



to obtain the mean {V2) 
standard deviation (Tdvn ---vt,-,, . 

elliptic flow parameter V2, where the factor comes from 
integration over 77. 



The induced V2 fluctuations arising from fluctuations 
in the number of participating nucleons are calculated by 
parameterizing the {V2) versus A^part results and folding 
them with the iVpait distributions in each centrality bin. 
The relative contribution of these fluctuations to crdyn 
is found to be less than 8%. Results in this letter are 
presented after subtraction of iVpart induced fluctuations. 

Systematic errors have been investigated in three 
main classes: variations to the event-by-event analy- 
sis, response of the analysis procedure to known input 
CTdyn, and intrinsic differences between HIJING events 
and data. Various modiflcations to the event-by-event 
analysis have been applied. Corrections, previously used 
in the hit-based event-plane analysis [51 [7], to account 
for signal dilution due to detector occupancy and to 
create an appropriately symmetric acceptance have been 
applied to both HIJING and data events. The thresholds 
for background hit rejection have been varied. These 
changes lead to at most 4% variations in the observed 
relative fluctuations demonstrating a good understand- 
ing of the response function. The determination of 
the response function and the flnal fitting procedure 
have been studied by performing the analysis on sets 
of modified HIJING events with varying input CTdyn- 
Differences between input and reconstructed (Tdyn are 
identified as a contribution to the systematic uncertainty. 
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FIG. 2: (Tdyn/{t'2) versus A^part for Au+Au collisions at 
— 200 GeV. The continuous and dashed thick black 
lines show (j{epart) / {(-part) calculated in Glauber MC [14] and 
CGC [15] models, respectively. The shaded grey band (for 
data) and thin black contour line (for Glauber MC) show 90% 
C.L. systematics errors. See text for discussion of comparing 
the plotted data to the models. 

The sensitivity of the measurement is observed to be 
Hmited for very low (U2) values. Therefore the 0-6% most 
central events, where the reconstructed (W2) is below 3%, 
have been omitted. Differences between HIJING and 
data in terms of dN/dr] and V2{ri) can, in principle, lead 
to a miscalculation of the response function. A sample 
of MC events has been generated, in which the dN/dj] 
distribution of HIJING events is widened by a simple 
scaling to match the measurements in data within the 
errors. The difference between results obtained with 
and without this modification, as well as the difference 
between results with two different parameterizations of 
V2 {f]) are identified as contributions to the systematic un- 
certainty. Other systematic studies include using a flat, 
rather than Gaussian, ansatz for the true V2 distribution, 
/(V2), and performing the analysis in different collision 
vertex and event-plane angle bins. The uncertainty in 
the contribution of A^part induced fluctuations has also 
been estimated via different parameterizations of the (^2) 
versus iVpart results. Contributions from all error sources 
described above are added in quadrature to derive the 
90% confidence level error. 

Fig. [1] shows the mean, (W2), and the standard devia- 
tion, Udyn, of the elliptic flow parameter V2 at midrapidity 
as a function of the number of participating nucleons, in 
Au-f Au colhsions at ^/s^ = 200 GeV for 6-45% most 
central events. The results for {V2) are in agreement 
with the previous PHOBOS V2 measurements |7], which 
were obtained with the event-plane method for charged 
hadrons within |?7| < 1. The uncertainties in dN/dr] and 
V2{rf), as well as differences between HIJING and the 
data in these quantities, introduce a large uncertainty 
in the overall scale in the event-by-event analysis due 



to the averaging procedure over the wide pseudorapidity 
range. The event-plane method used in the previous 
PHOBOS measurement s ar e known to be sensitive to 
the second moment, \/Jv^, of elliptic flow [l3]. The 
fluctuations presented in this letter would lead to a 
difference of approximately 10% between the mean, (W2), 
and the RMS, {^2), of elliptic flow at a fixed value of 
.^part- However, a detailed comparison is not possible 
for our (^2) measurements due to the scale errors, which 
dominate the systematic uncertainty on {V2) and ddyn- 
Most of the scale errors cancel in the ratio, adyn/{v2), 
shown in Fig. [2] revealing large relative fluctuations of 
approximately 40%. 

These results include contributions from both ellip- 
tic flow fluctuations and non-flow correlations. With 
no prior information on the direction of the reaction 
plane, it is not possible to disentangle these two con- 
tributions completely. However, several methods have 
been proposed to estimate the contribution of non- 
flow correlations to the observed dynamic V2 fluctua- 
tions. One can assume that the correlations in A-l-A 
collisions can be modeled by superimposing p-l-p col- 
lisions [22]. However, data from RHIC reveal many 
differences between the overall correlation structure in 
Au-|-Au and p-l-p (e.g. [^5H^ ). A more data-driven 
approach assumes that non-flow correlations will be small 
for particle pairs with large pseudorapidity separations 
(for example, A77 > 2) |57|. Under this latter assumption, 
it is estimated that the relative fluctuations in the actual 
elliptic flow account for a very large fraction (79-97%) 
of the observed relative dynamic fluctuations in the V2 
parameter [57] . No attempt was made to correct the 
data in Fig. [2] for non-flow effects since the validity of 
the large Arj assumption cannot be unambiguously tested 
with existing data. 

The measured dynamic fluctuations in V2 are directly 
comparable to models that incorporate both elliptic flow 
and two particle correlations. Furthermore, without 
making any assumptions about non-flow, these data 
establish an upper limit on the magnitude of under- 
lying elliptic flow fluctuations. Also shown in Fig. [2] 
are Cep^rt/i^part) at fixed values of Npart obtained in 
MC Glauber [l3] and color glass condensate(CGC) [TB] 
calculations. The 90% confidence level systematic errors 
for MC Glauber calculations (shown as a countour line 
in Fig. [2]) are estimated by varying Glauber parameters 
as discussed in Ref. [T^]. Due to the uncertainties in 
non-flow effects discussed previously, it is not possible to 
conclude which of these two models is more consistent 
with the measured dynamic V2 fluctuations. 

In summary, we have presented the first measurement 
of event-by-event V2 fluctuations in Au-f Au collisions at 
y/s^ = 200 GeV. The relative non-statistical fluctua- 
tions of the V2 parameter are found to be approximately 
40%. Independent estimates of the non-flow correlation 
magnitude suggest that the major contribution to these 
fluctuations are due to intrinsic elliptic flow fluctuations. 
We show that the magnitude and centrality dependence 
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of observed dynamic fluctuations are consistent with 

predictions for fluctuations of the initial shape of the 
collision region. These results provide qualitatively 
new information on the initial conditions of heavy ion 
collisions and the subsequent collective expansion of the 
system. 
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